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The possibility of driving an Anderson metal-insulator transition in the presence of scale-free dis- 
order by changing the correlation exponent is numerically investigated. We calculate the localization 
length for quasi-one-dimensional systems at fixed energy and fixed disorder strength using a stan- 
dard transfer matrix method. From a finite-size scaling analysis we extract the critical correlation 
exponent and the critical exponent characterizing the phase transition. 
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The Anderson model of localization^ has been subject 
of intense study over the past decades. In particular, the 
occurrence of a metal-insulator phase transition (MIT) in 
three dimensions (3D) has attracted a lot of interestPS 
Theoretical studies of the MIT have focused mainly on 
situations with uncorrelated disorder^HSl. Therefore, one 
of the open questions in the field is the role of long-range 
correlated disorder in the Anderson MIT. 

For uncorrelated disorder, the Anderson transition can 
be driven either by increasing the disorder strength or 
by changing the Fermi energy^. In the former case, for 
sufficiently strong disorder strength W > W C {E) all elec- 
tronic states are exponentially localized, where the value 
W C (E) depends on the Fermi energy E. On the other 
hand, at fixed disorder strength states with \E\ < E C (W) 
are extended and otherwise localized. 

The presence of correlations provides an additional 
possibility to achieve the MIT. Depending on the na- 
ture of the correlations, a transition may, in principle, 
be also driven by a change of the correlation strength 
or correlation length. Such a scenario might be relevant 
in situations where the disorder is induced by a complex 
environment surrounding the system of interest. 

In this Brief Report we study the possibility of a 
correlation-strength driven Anderson MIT in 3D. We 
consider the case of scale-free disorder, which is char- 
acterized by a power-law with correlation exponent a. 
We find at fixed energy and fixed disorder strength that 
the localization length behaves as 

X(a) oc \a c - a\~ v , (1) 

where the critical exponent v depends on the values of 
W and E. The obtained critical values a c are consistent 
with results for disorder and energy driven MITs in the 
presence of scale-free disorder^. 

To study the influence of scale-free disorder on the An- 
derson MIT, we use the usual tight-binding Hamiltonian 
in site representation^^ 

H = J>|i><i|-X>|i)<j|, (2) 

i ij 



where |i) denotes a localized state at lattice site i. The 
hopping matrix elements iy are restricted to nearest 
neighbors. As usual, we set these elements to one and 
thereby fix the unit of energy. The on-site potentials Si 
are taken as random numbers with a Gaussian probabil- 
ity distribution. Specifically, we use random potentials 
with mean (e{) = and a correlation function of the form 

C{i) = ( £i e i+£ ) cx \£\- a , (3) 

where a is the correlation exponent. In the con- 
text of Anderson localization, this correlation function 
has been used to study localization in the presence 
of long-range correlations for one-dimensionaPH51 two- 
dimensionaP^Sl anc j three-dimensionaF^Q systems. 

For the numerical calculations we generate the on-site 
potentials for systems of size M x M x L using a modi- 
fied Fourier filtering method (FFM|2H. Additionally, we 
shift and scale the resulting random numbers to have 
vanishing mean and variance C(0) = W 2 /12. We fo- 
cus on quasi-one-dimensional systems with L = 400000 
and M = 5,7,9,11 and 13. The localization length A 
is calculated using a standard transfer-matrix method 
(TMMp. Monitoring the variance of the change of the 
Lyapunov exponent during the TMM iterations gives a 
measure of the accuracy of the localization length^. We 
use a new seed for each parameter combination (E, W, 
a, M). Lastly, the critical exponent and the critical cor- 
relation strength are obtained from a finite-size scaling 
(FSS) analysis^. We expand the one-parameter scaling 
law for the reduced localization length A = X/M into a 
Taylor series 

A(M, r) = J2 ^M^F^xM 1 '") , (4) 

n=Q 

where x is a relevant scaling variable, <j) is an irrelevant 
scaling variable, y > is the irrelevant scaling expo- 
nent and r measures the distance from the critical point. 
However, instead of using energy or disorder strength to 
measure this distance, we utilize the correlation expo- 
nent, i.e., t = \a — a c \/a c . The functions F n , x and <f> 
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FIG. 1. Schematic phase diagram based on the results re- 
ported in RefP. The phase-space points discussed in the text 
are indicated by symbols (* and x). 



the band center. 

Qualitatively, the correlation-strength driven transi- 
tion can be understood by assuming an effective disor- 
der strength, W c s(a), which depends on the correlation 
exponent. An effective smoothening of the disorder po- 
tential has, for example, been observed for ID systems, 
where the localization length in the band center increases 
for smaller correlation exponent d 8 l 15 l It is also in ac- 
cordance with the shift of the phase boundary towards 
higher energies and stronger disorder shown in Fig.[l] Ac- 
cordingly, the transition occurs when Weff(a c ) = W C {E). 
Close to the transition, the localization length would di- 
verge according to A oc |Woff(a) — W c \~ v ° oc \a — a c \^" , 
where we have expanded the effective disorder strength to 
first order, W^f^a) « W c + (a — a c )dW c s/da\ ac . By con- 
struction this procedure yields the correct critical correla- 
tion strength, but it does not explain the deviation of the 
observed critical exponents from the universal value v . 
Provided the one-parameter scaling law holds in the pres- 
ence of long-range correlations, this discrepancy might 
also indicate that the FSS method in the normally used 



are further expanded up to order tir, ttir and mi, respec- 
tively. Taking n\ > allows us to consider corrections 
to scaling due to the finite size of the sample, which is 
reflected in a systematic shift of A with M in Eq. Q. 
Using a least squares fit of the expansion of the reduced 
localization length A to the numerical data allows us to 
obtain the critical parameters. Although one does ad hoc 
not expect the FSS analysis to be valid in the present 
case, we find that it is working surprisingly well, as we 
will show in the following. 

We set E = const, and W — const, while varying a. 
The chosen values of E and W are indicated in Fig. [T] 
which shows a schematic phase diagram for the Anderson 
MIT in the presence of scale-free disorder. The position 
relative to the transition boundaries provides a first esti- 
mate of the expected critical correlation exponent. 

In Fig. [2] the reduced localization length is shown in 
the vicinity of the band center, E = 0, setting W = 26. 
From the dependence on the system size M, a clear 
transition can be seen. For small correlation exponents 
(a < 1.5) the reduced localization length increases with 
increasing size M, while for large exponents (a > 1.5) 
it decreases. The former is characteristic for a metal- 
lic phase and the latter for an insulating phase. The 
FSS procedure yields for the critical correlation strength 
a c = 1.44 ± 0.04, which agrees very well with the value 
expected from the phase diagram. The critical exponent 
is v = 0.98 ± 0.09 (y = 2.0 ± 1.3), which is different 
from the value vo = 1.58±0.03 obtained for uncorrelated 
disorder^! and from via — 1.5) = 1.69 ± 0.22 reported 
for scale- free disorder^, both taken at E = 0. 

Also, for E = 6.0 and W — 16.5 we find a transition, 
as shown in Fig. [2j In this case the critical value is found 
to be a c = 2.85 ± 0.03, again consistent with the phase 
diagram in Fig. [T] The critical exponent is v — 0.62 ± 
0.03, which is even smaller than the exponent found at 
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FIG. 2. Reduced localization length A vs correlation exponent 
a. Solid lines show FSS fit to numerical data, (a) Taking 
corrections to scaling into account (tir = 2, nx = 1, itir = 
2, mi = 0) for E = 0.0, W = 26.0. (b) Without taking 
corrections to scaling into account (tir = 2, rii = 0, tur — 2, 
mi = 0) for E = 6.0, W = 16.5. 
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form is not suitable to extract the critical exponent in 
the present case. 

In summary, we have studied the influence of scale-free 
disorder on the Anderson MIT at fixed energy and fixed 
disorder strength. By varying the correlation exponent 



we found an increasing reduced localization length for 
a < a c and a decreasing reduced localization length for 
a > a c when increasing the system size. A FSS analysis 
yielded critical exponents which depend on the values of 
E and W and are smaller than the universal value isq 
found previously for uncorrelated disorder^!. 
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